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ABSTRACT 

In the theory of structure formation, galaxies are biased tracers of the underlying matter density 
field. The statistical relation between galaxy and matter density field is commonly referred as galaxy 
bias. In this paper, we test the linear bi as model with weak-l ensing and galaxy clustering measurements 
in the 2 square degrees COSMOS field (IScoville et al.l [2007). We estimate the bias of galaxies between 
rcdshifts z = 0.2 and 2 = 1, and over correlation scales between R = 0.2 h^ 1 Mpc and R = 15 hr 1 Mpc. 
We focus on three galaxy samples, selected in flux (simultaneous cuts Iguw < 26.5 and K s < 24), 
and in stellar-mass (10 9 < M* < 10 10 hr 2 M and 10 10 < M* < 10 11 hr 2 M Q ). At scales R > 
2 /j Mpc, our measurements support a model of bias increasing with redshift. The Tinker et al. 
(2010) fitting function provides a good fit to the data. We find the best fit mass of the galaxy halos 
to be log(M 20 o//i" 1 M Q ) = 11.7±°J and log(M2oo//i _1 MQ) = 12A±l l respectively for the low and 
high stellar-mass samples. In the halo model framework, bias is scaie-dependent with a change of 
slope at the transition scale between the one and the two halo terms. We detect a scale-dependence 
of bias with a turn-down at scale R — 2.3 ± 1.5h~ Mpc, in agreement with previous galaxy clustering 
studies. We find no significant amount of stochasticity, suggesting that a linear bias model is sufficient 
to describe our data. We use N-body simulations to quantify both the amount of cosmic variance and 
systematic errors in the measurement. 

Subject headings: cosmology: observations — gravitational lensing: weak — large-scale structure of 
Universe 



1. INTRODUCTION 

In the theory of structure formation, galaxies form in 
dark matter overdensities. However, the distribution of 
galaxies does not perfectly match the underlying matter 
distribution, and the relation between the two is called 
galaxy bias. 

The galaxy bias relation 5 g = f(S) relates the mat- 
ter density contrast S to the galaxy density contrast S g . 
The density contrast 5 is defined as the ratio between 
the local and the mean densities S = p/p — 1. To first 
order, the galaxy bias relation S g = bS is parametrized 
by the bias parameter b. However the stochasticity in 
the physical processes of galaxy formation might intro- 
duce some stochasticity in the galaxy bias relation. The 
amount of stochasticity can be measured with the cor- 
relation coefficient r defined as the ratio between the 
covariance and the variance of the density contrasts 5 

and S g , r — (S g 8)/ J (6 2 )(5 2 ). A correlation coefficient 

r = 1 means the bias relation is perfectly linear, whereas 
< r < 1 suggests stochasticity and/or non-linearity 
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(jDekel fc Lahavlll999D . We do not expect the bias re- 
lation to be anti-correlated with r < 0, nor b < (i.e. 
galaxy clustering decreases when matter clustering in- 
creases) . 

Bias varies for different galaxy populations. Galaxy bi- 
asing is known to be larger for luminous , red, and high- 
redsh i ft galaxies (iMarinoni et al.l 120051: [ Meneux et al.l 
120061 : iCoil et al.l 120081 iZehavi et all 12010 1 Bias also 
varies with scale, but differently for different galaxy types 
and luminosity (iCresswell fc Percivall 120091 : lOuchi et al.1 
20051 iMcCracken et alJl20TaT In the halo-model frame- 
work, the scale-dependence of bias is explained as the 
transition from a small-scale regime where galaxies 
are in the same halo of dark matter to a large-scale 
regime where gala xies are in two different halos (see e.g. 
iZheng et al.1 120091 for a recent discussion). The scale- 
dependence of bias therefore tells us about the dark mat- 
t er halo properties arou n d gala xies . 

iTegmark fc Bromlevl (|1999f ) first measured bias 
stochasticity in the Las Campanas Redshift Survey. At 
the same time, simulations started to predict stochastic- 
ity, as well as scale and redshift-de pendence in the bias 
relation due to physical processes (IBlanton et al.1 [l999t 
iSomerville et alj I2001t lYoshikawa et al.l I2001D In the 



2dF Galaxy Redshift Survey. IWild et al.1 (12005;, used 



the count-in-cell technique to measure the non-linearity 
and stochasticity of the bias relation for red/early 
and blue/late- type galaxies at redshift z < 0.114. For 
both categories, they found more stochasticity at scale 
R = 10 Mpc than at scale R = 45 Mpc. Recent results 
from simulations also predict an increasing amou nt of 
stochasticity at smaller scales (jBaldauf et al.l[2010f ) . 
The correlation coefficient is at the same time an esti- 
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mate of stochasticity a nd non-linearity. I n the zCOS- 
MOS redshift survey (jLillv et all 120091 ) . iKovac et all 
(2009]) found non-linearity to contribute less than ~ 
0.2% to the bias relation for luminosity-selected galax- 
ies (Mb < —20 — z). Therefore in the following, we will 
assume that the correlation coefficient r mostly quanti- 
fies bi as stochastici t y. W ith respect to the bias param- 
eter b, IKovac et al.l (|2009D found the bias relation to be 
scale-independent between scales R = 8 hT x Mpc and 
R = 10 hr 1 Mpc, and redshift-dependent between red- 
shifts z — 0.4 and z = 1. 

In this wor k, we measure ga laxy bias with a technique 
proposed by iSchneideil (|1998| ) based on galaxy cluster- 
ing and weak lensing. In short, weak lensing is used to 
derive matter clustering from the shape of background 
galaxies, which is then compared to the gal a xy clu ster- 
ing to infer the bias. Ivan Waerbeke et al.l (jl999H em- 
phasized the potential of this technique with analyti- 
cal calcu lations for surveys with vario us depths and ar- 
eas, and iHoekstra et~aH (|200lL 12002(1 applied it to the 
50.5 deg 2 of the RCS and VIRMOS-DESCART sur- 
veys. They measured the linear bias parameter b and 
the correlation coefficient r for a flux-limited galaxy sam- 
ple (19.5 < Rc < 21) at redshift z ~ 0.35, and on 
scales between R = 0.2 and 9.3 Mpc. They found 
strong evidence that both b and r change with scale, and 
r ~ 0.57 at 1 h§Q Mpc, which suggests a significant de- 
gree of stoch asticity and /or non -linearity in their sample 
of galaxies. iSimon et al.l (|2007f ) applied the same tech- 
nique to the 15 deg 2 of the GaBoDS survey for three 
flux-limited galaxy samples in the R-band, at redshift 
z =0.35, 0.47 and 0.61. They also found bias to be scale 
dependent, with an increasing amount of bias at small 
scales. In addition, they found the bias parameters b 
and r to be redshift-independe nt within statistical u ncer- 
tainties, with r ~ 0.6. Finally. TSheldon et al.l (j2004H per- 
formed a comparable study in the 3800 deg 2 of the SDSS 
survey, with volume- and magnitude-selected galaxies 
(0.1 < z < 0.174 and -23.0 < M r - 5 l og h < -21.5). In 
agreement with Hoekstra et al. (20Q2|), they found the 
ratio b/r to be scale-independent between 0.4 and 6.7 
h- 1 Mpc, with (b/r) = (1.3 ± 0.2 )(n m /0.27) for a flux - 
limited galaxy sample. In fact, the lHoekstra et all (|2002[ ) 
results suggested that the respective scale-dependence of 
b and r conspire to produce a scale-independent ratio 
b/r. In summary, the consistent value of r < 1 found in 
all these studies suggests some degree of stochasticity in 
the investigated galaxy samples. 

In this work, we study the bias of galaxies in the 2 
square degrees of the COSMOS field. We test the lin- 
ear bias model, and its dependence on scale and redshift. 
For two galaxy samples selected in stellar mass and one 
sample selected in flux, we measure bias on scales be- 
tween 0.2 /i" 1 Mpc and 15 h" 1 Mpc, and redshifts be- 
tween z = 0.2 and 2 = 1. Thus, we extend the pre- 
vious works on bias both in scale, redshift and stellar 
mass. COSMOS is ideal to perform this study because 
accurate stellar mass have b een measured up to redshift 
z ~ 1.2 (jBundv et al.l |2010|) . and photom etric redshifts 
up to redshift z ~ 4 (jllbert et aLl|2009f l. In addition, 
with the high resolution images from the the Advanced 
Camera for Survey aboard the Hubble Space Telescope 
(ACS/HST), we can precisely measure the shape of more 



than 200,000 galaxies. Measuring the shape of galaxies 
at redshift z > 2 is crucial for weak lensing to allow ac- 
curate determination of the matter density up to redshift 
z ~ 1. In section [21 we present the COSMOS survey, our 
lensing catalog and our selection of foreground galaxies 
for which we want to measure the bias. In section [3J we 
describe the lensing formalism used in this paper. In sec- 
tion |U we present a set of simulations we use to quantify 
the impact of cosmic variance. In section [5j we show our 
results, and we go on to discuss some possible systematic 

errors in section [6] We conclude in section [7] 

In this paper, we assume the WMAP7 (jKomatsu et al.l 
2010*) cosmological parameters f2 m = 0.272, Q\ = 0.728, 
fl b = 0.0449, cr 8 = 0.809, n = 0.963, u = -1 and 
h = 0.71. 

2. DATA 

2.1. General properties 

The COSMOS field is a 1.64 deg 2 patch of sky close 
to the equator (a=10 h 00 m 28?, S=+2°l2'.21 'f). It was im- 
aged with ACS/HST during c ycles 12-13 (jScoville et al.l 
2007; Ko ekemoer et al.l 12007( 1. An intensive follow-up 
from UV to IR was conducted by several ground-based 
and space-based telescopes. In this study, we take ad- 
vantage of the spatial resolution of the ACS imaging 
to measure the shape of the galaxies, and of the multi- 
filter ground- based imaging to d erive accurate photomet- 
ric redshifts (jllbert et all 12009( 1. 

2.2. Lensing catalog 

We use an updated version of th e lensing catalog 
presented in ILeauthaud et al.l (|2007l ). In this new 
catalog, raw ACS/f814w images have been corrected 
for the Charge Transfer In efficiency (CTI ) effec t using 
the algorithm described in Ma ssev et al.l (j2010H . Ob- 
jects are d etected with the "hot-c old" technique de- 
scribed in ILeauthaud et all (|2007[ 1. i.e. Sextractor 
(jBertin fc Arnoutslll996( l is run twice : first to detect ex- 
tended objects, and then compact objects. Regions con- 
taminated by bright stars or satellite trails are masked 
out. Point sour ces and spurious detect ions are discarded 
as detailed in (jLeauthaud et al] l2007f) . The catalog is 
complete up to magnitude MAG_AUTO = 26.5 in the 
ACS/f814w band (hereafter Jguw)- 

Galaxy ellipticiti es are measur e d with the RRG algo- 
rithm described in iRhodes et al.l (|2000[ K and calibrated 

in [L 



with t he same shear calibration as in iLeauthaud et al] 
(20071). In 

our final source catalog, we make the cuts on 
size d > 1.6 times the PSF size and S/N > 4.5 as in 
ILeauthaud et al.l ((2007). We also remove galaxies with 
more than one peak in their photometric redshift prob- 
ability distribution function, or galaxies masked in the 
Subaru imaging. This way, we obtain an average redshift 
error er z /(l + z) < 0.05 Our source catalog with good 
shapes and redshifts contains 210,477 galaxies, which 
yields a density of 36 gala xies per square arc-minute. 
As in Ma ssev et al.l (|2007bj l. we assign to every source 
galaxy a inverse variance weight w aimed at maximiz- 
ing the S/N, and calculated from the galaxy apparent 
magnitude 



1 



w 



0.1 ' 



(1) 
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where the observational noise (pixellation, shot noise, 
etc) is well approximated by the following function of 
the galaxy apparent magnitude in the Isuw band 



a e = 0.32 + 0.0014(MAG_AUTO - 20) 3 



(2) 



The second term 0.1 represents the variance of the intrin- 
sic shape noise, which w as found to be pretty constant in 
ILeauthaud et al.l (120071). Note that the p ower of 2 was 
omitted as a typo in lMassev et~aLl (I2007bfl . 

2.3. Foreground and background galaxy catalogs 
2.3.1. Catalog cuts by redshift 

In this work, we want to study the evolution of galaxy 
bias as a function of redshift and scale. Bias is computed 
from the ratio of galaxy and matter clustering in bins of 
redshift. To ease the interpretation of the bias results 
later on, we choose to define foreground galaxy samples 
with redshift distributions matching the redshift distri- 
bution of the matter that most efficiently perturbs the 
background galaxy shapes through weak lensing. Typ- 
ically, the matter that lies halfway between us and the 
background galaxies, is the one that most efficiently per- 
turbs their shape. We show in Section [5] that a mismatch 
in redshift can change the estimated bias. 

We create 3 bins Bl, B2 and B3 of background sources. 
We set the bin limits to z = 0.3, 0.8, 1.4 and 4, so that 
the number of galaxies and signal to noise in each bin 
is similar. We assign a galaxy with redshift z ± cr 6 s% 
to a particular redshift bin with limits z\ ow and Zhi g h if 
z — °~68% > z \ow and z + er 68 % < z h i g h- The values of 
z and <J68% are estimated from the photometric redshift 
likelihoods. The number of galaxies per bin is reported 
in Tabled! 

For each background bin, we compute the lensing ef- 
ficiency. The lensing efficiency is a calculation based on 
the redshift distribution of the background galaxies in co- 
moving coordinates Pb(w) = pb(z)dz / dw and of the cos- 
mological parameters. We use it as a weighting function 
in the following to project the 3D matter power spec- 
trum along the line of sight. The lensing efficiency as a 
function of the comoving distance to the lens is defined 
as 

2c z a(w) J w f K w' 

where the functions a(w) and fi{(w) are respectively the 
comoving scale factor and the comoving angular distance. 
We find the lensing efficiency curves to peak at redshifts 
Zbi = 0.22, z B 2 = 0.37 and z B3 = 0.51. 

Next, we create 5 bins of foreground galaxies. For 
the three first bins, we adjust their centers to match the 
peaks of the 3 lensing efficiencies. Then, we adjust their 
width so that the signal to noise of the galaxy-galaxy 
cross-correlation measurements is the largest. Too broad 
bins increase the mixing of angular and physical scales, 
and too narrow bins increase shot noise. Bins must 
also be broad enough, so that we can use the Limber 
approximation to compute the theoretical signal. In- 
deed, the Limber approximation br eaks down at large 
scales for too narrow redshift bins. ISimonl ((20071) com- 
puted the minimal bin width as a function of scale and 



redshift, beyond which the Limber approximation be- 
comes inaccurate by more than 10%. To minimize the 
impact of this issue, we fix the bin width so that we 
reach 10% inaccuracy at 100 arc- minutes, i.e. twice the 
largest scale at which we compute the correlation func- 
tions. In summary, we obtain three foreground bins with 
limits 0.17 < z F1 < 0.27, 0.32 < z F2 < 0.42, and 
0.47 < zp 3 < 0.57. To these 3 foreground bins, we 
add two extra bins to probe the bias at higher redshift 
0.60 < z Fi < 0.75 and 0.87 < z F5 < 1.07. Finally, we as- 
sign the foreground galaxies to their respective bins using 
the same criteria we used for the background sources, and 
compute Pf(z) the redshift distribution of the foreground 
bins. Again, we will use these redshift distributions as 
weighting functions to project the 3D matter power spec- 
trum along the line of sight, and compute the different 
estimators presented below. 

Note here three important points regarding our anal- 
ysis. First, we compute the redshift distribution of the 
galaxies that fall in the foreground and background red- 
shift bins Pf{z) and Pb(z) by summing their redshift pos- 
terior probability distribution function (pdf) p( % ' (z). The 
redshift distributions Pf(z) and Pb(z), as well as the lens- 
ing efficiency curves are shown in Figure [TJ Summing the 
redshift pdfs instead of building histograms of the indi- 
vidual redshift estimates has several advantages : (i) it 
helps to deal with skewed pdfs, (ii) it tells us about the 
probability of having galaxies outside a redshift bin. In 
spite of our cuts in redshift, we still find that 30% of 
the galaxies are probably outside the redshift limits of 
their respective redshift bins. If instead we choose to cut 
galaxies at 3<768% from the redshift bin limits, we only 
decrease this probability to 28%. (iii) It gives us theoreti- 
cal predictions that better reproduce the measurements. 
Summing the pdfs produces redshift distributions with 
tails falling outside the redshift bins limits. As a result, 
we predict projected power spectra with larger ampli- 
tudes at small and large scales. They reproduce better 
the measurements. 

Second, we must emphasize that the redshift distribu- 
tions we obtain are weighting functions that we use to 
project the 3D matter power spectrum along the line of 
sight, and predict the behavior of several estimators. In 
particular, we use the lensing efficiency in place of the 
lensing-effective matter distribution. This latter would 
be the true matter redshift distribution multiplied by the 
lensing efficiency. Nonetheless, by averaging over many 
line of sights, we expect the lensing efficiency to approx- 
imate the lensing-effective matter distribution, and thus 
produce accurate predictions. 

Third, the centers of foreground bins F4 and F5 do 
not perfectly match the peak of the lensing efficiency for 
catalog B3. In section [BJ we show how such a mismatch 
can alter the bias measurements. We show that the ef- 
fects are particularly severe at small scales in the non- 
linear regime. Fortunately for bins F4 and F5 most of 
the angular scales we probe correspond to physical scales 
in the quasi-linear and linear regime. Since the scale- 
dependence of the bias in the linear regime is small, the 
effect of mismatch in redshift will be small. 



2.3.2. Catalog cuts by flux and stellar mass 
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TABLE 1 

Redshift range and number of galaxies per foreground 
and background bins. 



Bin Redshift 



Stellar Mass Range [/i 72 2 M Q j 





Range 


10 9 -> 10 10 


10 10 -> 10 11 


All 


Fl 


[0.17, 0.27J 


1,298 


644 


7,920 


F2 


[0.32, 0.42] 


3,151 


1,770 


13,898 


F3 


[0.47, 0.57] 


3,215 


1,456 


11,256 


F4 


[0.60, 0.75] 


7,776 


3,899 


22,958 


F5 


[0.87, 1.07] 


11,828 


6,223 


27,911 


Bl 


[0.3, 0.8] 






63,529 


B2 


[0.8, 1.4] 






70,157 


B3 


[1.4, 4.0] 






59,056 



In order to investigate the dependence of bias with stel- 
lar mass, we further partition our foreground catalogs 
into one flux-selected and two stellar-mass selected cata- 
logs. The COSMOS catalog for the foreground galax- 
ies is complete u p to I^4\v < 26. 5 and K K < 24 

in AB magnitude (iCapak et al.l l2007t lllbert et all 120091 : 
iLeauthaud et al.ll2007l )~ The limit in X s -band is more 
conservative and comes from estimating the observed 
magnitude of a maximal M*/L stellar population model 
with solar metallicity, no dust, and a r = 0.5 Gyr burst 
of star fo rmation at re d shift Zf nr m = 5 (jLeauthaud et all 
l2011bH . iBundv et all <f2010h computed the 80% com- 
pleteness in stellar mass for magnitude selected galaxies 
K s < 24 and Isuw < 26.5. He found the 80% complete- 
ness lower limits in log M*/Mq to be 8.8, 9.3, 9.7 and 
10.0 for his redshift bins [0.2,0.4], [0.4,0.6], [0.6,0.8] and 
[0.8,1.0]. In order to have enough galaxy per bin, we bin 
in the two stellar mass bins 10 9 < M* < 10 10 h? 2 M 

and 10 10 < M* < 10 11 frf 2 2 M . Our low stellar-mass 
galaxy samples is therefore more than 80% complete in 
all redshift bins, and our high stellar-mass galaxy sam- 
ple is more than 80% complete in bins Fl, F2, about 
80% complete in bins F3, and less than 80% complete in 
bins F4 and F5. For comparison, the mean stellar mass 
in our flux-limited galaxy samples in bins Fl, F2, F3, 
F4 and F5 are 9.6, 9.7, 9.8, 10.0 and 10.1 in units of 



log(M*//if 2 2 M Q ). 

3. LENSING METHOD 

In gravitational lensing theory, the mass of foreground 
structures locally perturbs space-time, deviates the path 
of light-rays coming to us from background galaxies, and 
distorts their intrinsic shape. By analyzing the shape 
of many background galaxies, it is therefore possible to 
infer the properties of the intervening structures. For 
instance, the correlation of ellipticity measurements of 
background galaxies tells us about the clustering of the 
foreground structures. As such, gravitational lensing is a 
direct probe of the matter density field projected along 
the line of sight. Perturbations of the background galaxy 
shapes are of two kinds. The convergence k quanti- 
fies the amount of isotropic amplification, and the shear, 
7 = 71 + i^i quantifies the amount of stretch along some 
direction. 

If we consider a sample of background galaxies with 
lensing efficiency as a function of redshift g(z), the con- 
vergence and shear produced by the matter along the line 
of sight, expressed in terms of density contrast S is 



k(9) 



dw g(w)6(f K {w)9,w) . 



7(0) 



(4) 



(5) 



iSchneider et al.l (|l998) introduce the mass aperture 



statistic M ap defined as 



M ap (9) = / d 2 dU(\d\) K (d), 
Jo 



(6) 



where the integral extends over a disk of angular radius 9. 
U(9) is a compensated weight function which vanishes for 
i9 > 9. In a similar manner, the aperture number count 
of galaxies at a scale 9 is defined as 



N(6)= / d 2 4U{\4\)8 g (4) 
Jo 



(7) 



where U{d) is the same function as above, and S g (9) is 
the galaxy density contrast. 

3.1. Variance aperture statistics 

ISchneider! (|1998l) and Ivan Waerbeki (|1998D developed 
a formalism based on aperture statistics to measure the 
dependence of bias with scale. They define the follow- 
ing aperture statistics from the auto- and cross-power 
spectra of galaxies and matter P n {£), P K (£) and P nK (£) 
respectively 



(N 2 (9))=2n di I P n (l)[I(i9)Y 



(ML(6)) = 2n deep K (e)[i(£9)Y 



(A/"(0)M ap (0)) = 2tt / dlt P nK (£) [I{W)f 



(8) 

(9) 
(10) 



The filter I(x) is defined as (ISchneider et all 119981 : 
iHoekstra et~aTll2002l: [Simon et al.ll2007ft 



I(x) 



12 Ji{x) 



(11) 



where Ji(x) is the Bessel function of order 4. This filter 
is a narrow band filter, which peaks at x ~ 4.25. The 
power-spectra P n (£), P K {(-) and P Kn (£) derive from the 
3D matter power-spectrum P m {k, w) projected along the 
line of sight, weighted by the functions g(w), Pf{w) or 
Pf(w)g(w), and altered by the bias parameters b and r. 
Their respective analytical expressions are 



Pn{i) - 
P K {i) 

and 

Pkti (^) — 



L 



.. dW [f K (wW Pm 



Uh dw 9 ^ )Pf{w J brP„ 



fK{w)' 



Ik{w) 



,w 



(12) 
(13) 

(14) 
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Redshift 

Fig. 1. — Redshift distributions of the 3 background catalogs Bl— B3 (purple/dark line), and the five foreground catalogs F1-F5 (yel- 
low/light line). The normalized lensing efficiency curves (red dashed lines) give the location where the matter most efficiently perturbs the 
shape of the background galaxies. [See the electronic edition for a color version of this figure.] 
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Fig. 2. — Bias calibration factors fl and fl used to calibrate 
b and r respectively. We compute them assuming a lSmith et al.1 
(2003) non-linear p ower-spectrum w i th unb iased foreground galax- 
ies (6 = r = 1), a Eiscnstcin & Hu (f999) transfer function , and 
the WMAP7 cosmological parameters (Komatsu ct al. 2010). De- 
viations from a straight line are due to mismatches in Fourier space 
of scales probed by the galaxy and the matter aperture statistics 
(N 2 (8)) and (M^ p (9)) respectively. A straight line would indicate 
a perfect match between the two statistics. 

We use the Limber approximation here to simplify the 
integrals. 

The technique developed in this paper is powerful be- 
cause all the power spectra are filtered by the same nar- 
row band filter I(x). In contrast, the estimators £+, 
us (9) and (74) involve four different filters (see Ap- 
pendix IB]) . Combining them directly would lead to a 
mixing of modes in Fourier space, and would hamper 
the interpretation. 

3.2. Bias measured from aperture statistics 

The bias parameters b and r are defined as ratios of 
the variance aperture statistics 



&(0) = /i(0,O m ,f2 A ) x 



(Mi p (9)) 



r{9) = f 2 (e,n m ,Q A ) x 



(Af(9)M ap (6)) 



(15) 



(16) 



The functions f% and fi correct for the fact that dif- 
ferent cosmological volumes are probed by the different 
statistics. They are defined as 



W) = 



(AP(0)) 



=b=l 



h{6) 



(Ni{6)){MU6)) 



{N{9)M ap {6)) 



(17) 



(18) 



r=b=l 



and are computed assuming a ISmith et all |2003) non- 
linear power-spectrum with unbiased for e ground galax- 
ies (b = r = 1), a lEisenstein fc Hul (|1999|) trans- 
fer function, and the WMAP7 cosmological parameters 
(jKomatsu et alJl2010t ). 

In Figure O we show the evolution of fi and fi as 
a function of scale and redshift. These functions are 
not constant in scale. The amplitude of the deviations 
is larger for the higher redshift samples of galaxies. In 
Figure [21 function fi for pair Fl-Bl is nearly constant, 
whereas it varies by 25% for pair F5-B3. The scale- 
dependence of fi and fi is due to the fact that many dif- 
ferent physical scales are projected into the same angular 
bins, and the weighting functions in {N 2 (9)), (M 2 p (9)) 
a nd (N{9)M»J9)) are different . 

Ivan Waerbekd (|1998l ) showed that f\ and fi do not de- 
pend on the shape of the power-spectrum. We performed 
a set of tests to verify this statement, and reached the 
same conclusion. We compute d fj and fi with a power 
spectrum fo r the linear regime iPeacock fc Doddsl (|1996| ) 
instead of a ISmith et al.l ([2003D power-spectrum. We al- 
tered the amplitude cr 8 and spectral index n s by 20%, and 
we measured less than a 5% difference. The amplitude of 
the power spectrum does not have a strong influence on 
fi and f%. Note howeve r that the calibrat ion functions 
fi and fi scale with fl m (jSimon et al.ll2007l) . We fixed it 
to the WMAP7 value. 

3.3. Physical scales 

Because we make measurements in bins of angular sep- 
arations rather than physical separation, the bias param- 
eters we measure are averages of the true bias parameters 
b(k,w) and r(k,w) in k-space and com oving distances. 
Their respective expressions are given by ( Hockstra et al. 
12001) 

C h dw hi(8, w)b 2 (■qt^ttt , w) 

W {9) ~ fmT, 7—r Q \ > ( 19 ) 

J aw hi(a, w) 
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-)(<?) = 



Caw h3 (e,w)r{j^,w) 

f™ h dw h 3 (0,w) 



(20) 



where the weighting functions hi(0,w) and h 3 (9,w) are 
respectively 



hi(e,w) 



Ik{w) 



(21) 



ftaCO, «0 = P ^ W } 9 [^ PmUe, w) , (22) 
and the filtered matter power spectrum is 



Pfilter(0, W)=2K d££ P, 



,W [!(&)] 



fn{w) 

(23) 

Note that in Eq . [T9l and l20l the bias parameters b(k, w) 
and r(k,w) are estimated at scale k = 4.25/#/k-(io). In 
these equations, it was assumed that the bias parame- 
ters are evolving slowly with scale, and their product by 
the I(x) filter was equivalent to the product by a Dirac 
8d{x) function. Note as well that the functions Pf{z) 
and pf{z)g{z) are narrow in redshift space. As a result, 
very few angular scales are mixed together when bias pa- 
rameters b(k, w) and r{k, w) are projected along the line 
of sight. Therefore, the shapes of the bias parameters (b) 
and (r) in real space and b(k, w) and r(k, w) in k-space 
are very similar. 

We use the weighting functions hi and h 3 to relate 
the physical scales to the angular scales at which bias 
parameters are measured. The scales in physical units 
are obtained via the following expressions 



(R b (0)) = 2tt 



(R r (0)) = 2tt 



f™ h dw h\(9, w) 
Cdwh^Mef^; 

C dwh 3 (8,w) 
Cdwh 3 (9,w)ef^ 



(24) 



(25) 



In Figure |3l we show the weighting functions h\ and h 3 
for our 5 different pairs of redshift bins. The bias we mea- 
sure in the following is the averaged bias of the galaxies 
in the hatched regions. The h± and h 3 functions extends 
over the whole redshift range, but we only hatched the 
regions containing 68% of integrated signal in h\ and h 3 . 
We found the widths of the hatched regions to be scale 
independent. 

4. SIMULATIONS 

We use N-body simulations to quantify the effect of 
cosmic variance on our bias measurements. Theoretical 
estimates of the cosmic variance for mass aperture statis - 
tics exist ([Schneider et all l2002t Uoachimi et~aTl |2008), 
but they break down in the non-linear regime, making 
simulations a necessity. These simulations will also help 
us to check our tools. Indeed, the mock galaxies in these 
simulations are pure dark matter particles, and as such 
perfect tracers of the underlying dark matter. With our 
tools, we should therefore find b = r = 1 at all scales, 




Fl-Bl F2-B2 F3-B3 F4-B3 



□ 



F5-B3 Z) 

m 



h 3 (0, z) 



mi 



0.0 



0.4 0.8 
Redshift 



1.2 



Fig. 3. — The 5 weighting functions h\(9,z) and 113(8, z) used 
to compute the analytical predictions of the bias parameters (6) 
and (r) for each bin Fi-Bj. In both panels, the 5 hatched areas 
contain 68% of the integrated signal in h\ and /13. The width of 
the hatched regions is scale-independent. 
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Fig. 4. — The error budgets for the galaxy bias parameters (b)(8) 
and (r)(0), obtained for foreground Fl and background Bl. Shot 
noise and cosmic variance are obtained from simulations, whereas 
shape noise comes from the real data. 

and all redshifts. Any deviation would highlight a flaw 
in our tools, or an artifact in the method. 

We generate a set of 7 light-cones of 10x10 square de- 
grees each, and extending from redshift z — to z — 4. 
The particle mass is 7.5 x 10 10 M©. We estimate the 
shear and convergence on a 3D grid (fine in the angu- 
lar direction, coarse in the radial direction), using a no- 
radial-binning method (Kicssling ct al. 2010) . 

We make sure that the redshift distrib ution of the par- 
ticles matches a lEfstathiou et al.l ()1991h redshift distri- 
bution p(z) oc z 2 exp[— (z/zq) 2 ], with zq = 0.78 and 66 
galaxies per square arc-minute. This redshift distribu- 
tion is close to the actual redshift distribution in COS- 
MOS (jMassev et al.ll2007bl ). Finally, each simulated cat- 
alog is cut into 36 COSMOS-size catalogs, and we remove 
masked galaxies in the same way as we do with the real 
data. We obtain 343 mock catalogs of about 400,000 
galaxies each. We bin these galaxies in redshift, and ran- 
domly resample them so that the number of foreground 
and background galaxies match the numbers in Table [TJ 
The redshift of the galaxies is computed analytically from 
their comoving distance assuming a WMAP7 cosmology. 
To keep it simple, we do not simulate photometric red- 
shift uncertainties. 
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Figure 0] shows the error budget for the bias b and 
the correlation coefficient r computed for the foreground 
rcdshift bin Fl, and the background bin Bl. We split the 
error budget into cosmic variance, shot noise and shape 
noise. To quantify the uncertainty in our results due 
to cosmic variance, we compute the bias parameters for 
each COSMOS realization. The variance of the resulting 
distribution is therefore due to cosmic variance and shot 
noise. To quantify the amount of shot noise only, we 
bootstrap the galaxies in one single COSMOS realization. 

As shown in Figure EJ we find the cosmic variance for 
parameter (b) computed in foreground bin Fl and back- 
ground bin Bl to be larger than the shape noise. In con- 
trast for the correlation coefficient (r) , shape noise dom- 
inates at all scales, (r) is more sensitive to shape noise 
because it contains the contributions from (Af(6)M ap (6)) 
and (Mf (#)), which are both affected by shape noise. 
We perform the same analysis with the other pairs in 
bins Fi-Bj, and find that shape noise is always larger 
than cosmic variance both for bias parameters (b) and 
(r). The reason is that background galaxies in bins B2 
and B3 are fainter, so their shape estimation is more 
noisy. 

We also used the simulations to test our treatment 
of the photometric redshift, and especially the way we 
reconstruct the weighting functions pt(z) and Pf(z) for 
the background and foreground redshift bins. We sim- 
ulated photometric redshifts with different redshift un- 
certainties from d z /(l + z) = 0.01 to d z /(l + z) = 0.05. 
We reconstructed the weighting functions of bins Fl and 
Bl. For bin Fl with d z /(l + z) = 0.05, we measured 
the rms between the recovered and the true p(z) to be 
rms = 4.49 ± 0.32 for the standard histogram technique, 
and rms = 4.10 ± 0.33 for the sum of pdfs technique. 
For bin Bl with errors d z /(l + z) = 0.01, we mea- 
sured rms — 0.48 ± 0.01 and rms = 0.65 ± 0.01 re- 
spectively. Therefore, our technique is better for large 
rcdshift uncertainties, and lower redshift bins. In our 
particular case, where redshift uncertainties are about 
d z /(l + z) — 0.05 for the faint part of the samples, our 
technique is globally as good as constructing a standard 
histogram of the best fit estimates. As already mentioned 
previously, we preferred this technique because galax- 
ies outside the redshift limits were properly taken into 
account, which resulted in a better agreement between 
predicted and observed signals. To push this compari- 
son further, more realistic galaxy redshift pdfs should be 
simulated, but this is out of the scope of this paper. 

Note finally that our estimates of the systematic er- 
rors in our measurements of (6) and (r) might be un- 
derestimated because we do not consider scatter due to 
inaccurate photometric redshifts, and we do not popu- 
late the dark matter simulations with a realistic model 
of galaxies. Therefore, any additional scatter introduced 
by physical processes involved in galaxy formation will 
not be taken into account. 

5. RESULTS 

In this section, we present our measurements of the cor- 
relation functions. We use these measurements to derive 
the bias, and estimate its scale- and redshift-dependence. 

5.1. Correlation functions 



As described in Appendix [A] aperture statistics are 
computed from auto- and cross-correlation functions. We 
use the software Athen^H, which implements a tree- 
code, and computes auto- and cross-correlation func- 
tions. For the tree-code, we choose an opening anglc0 
of 0.04 radians. All correlation functions are measured 
in 950 logarithmic bins between 0.05 and 50 arc-minutes, 
which corresponds to a sampling scale of Az9 = 0.007i9. 
The aperture statistics and the bias parameters are also 
computed with this fine binning. We estimate the error 
bars with bootstrapping, i.e. we repeat the measure- 
ments 200 times with different galaxy samples randomly 
drawn from the input catalogs. For the figures, we rebin 
the data in coarse bins using a median binning technique. 
We used the median rather than the mean because the 
distributions of the bias parameters are not Gaussian, 
and so the median is more robust. However, this is a 
small correction since with the mean instead, the bias 
parameters b and r would be only about 3% larger. To 
allow a straightforward comparison, we use the same bin- 
ning for both the data and the simulations. 

In the following, we will always co mpare our results to 
numerical predictions. We found the S mith et al.l (|2003T ) 
fitting formula to give a good approximation of the mea- 
surements at large scales. In contrast at small scales, we 
found the fitting formula to under- predict the measure - 
ments (see e.g. Figures l6l [8l and f9]> . Iffllbert et al.l (2009) 
had already noticed this discrepancy at scales i > 10000, 
and suggested that it could be due to th e low r esolution 
of the simulations used by ISmith et al.l ([2003D . In our 
case, the discrepancies could at least partly also originate 
from the fact that we compare power-spectra of biased 
galaxy samples and dark-matter. Indeed, power-spectra 
of galaxies are expected to deviate from pure dark mat- 
ter power-spectra. Analyzing these deviations as a func- 
tion of scale, redshift and galaxy properties is in fact the 
whole purpose of this paper. Therefore, comparing the 
measured and the predicted measurements already give 
us a first guess by eye of the amplitude of galaxy bias for 
a particular galaxy population. 

Galaxy clustering — In Figure [5] we show the pro- 
jected galaxy auto-correlation function uj(8) for our 5 
foreground redshift bi ns. We compute uj(6) using the 
lLandv k, Szalavi (|1993| ) estimator defined as 



, ns DD DR 

oj(0) = — 2 

RR RR 



1 



(26) 



where the normalized numbers of pairs DD, RR and DR 
refer to pairs of galaxy positions, random positions, and 
galaxy and random positions respectively. 

By bootstrapping the data, we include shot noise in 
the error bars for the data, and cosmic variance for the 
simulations. At high redshift, we calculate smaller error 
bars because (i) high redshift bins contain more galaxies, 
and (ii) we probe larger volumes and cosmic variance de- 
creases in larger volumes. In bins F4 and F5, we observe 
2cr and 3cr deviations respectively between the measured 

7 http:/ /www2. iap.fr/users/kilbinge/athena 

8 In Athena, galaxy properties are spatially averaged in nodes. 
The opening angle between a central and a tangential node is the 
ratio between the size of the tangential node divided by the distance 
between the nodes. The properties of two nodes are correlated if 
the opening angle is smaller than 0.04 rad. 
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Fig. 5. — Projected galaxy auto-correlation for our flux-selected 
galaxy samples in bins Fl to F5. (dashed line) Fit to the data 
with a power-law. The best fit scaling parameter A w and slope 
8 are reported in Table [2] (red line) Power-law with parameters 
inferred from the fit to (Af 2 (0)). This represents the signal we 
should obtain after IC correc tion of u>(9) (see te xt), (bold line) 
Predicted signal derived from a lSmith et al . ( 2003 ) power-spectrum 
with unbiased galaxies. Error bars include shot noise (thick line) 
and cosmic variance (thin line) . Cosmic- variance always dominates 
the measurements. 



TABLE 2 

Scaling and slope of the power-law fit to u>(9) 





Direct fit 


IC Corrected 


using (M' z (e)) 


Bin 


A W {V) 


S 


A W (V) 


5 


IT 


0.32 


0.75 


0.61±0.28 


0.37±0.12 


F2 


0.21 


0.93 


0.50±0.16 


0.41±0.09 


F3 


0.15 


0.93 


0.22±0.05 


0.64±0.08 


Fl 


0.18 


0.65 


0.21±0.03 


0.54±0.05 


F5 


0.14 


0.63 


0.15±0.02 


0.55±0.03 



Note. — The IC corrected amplitu des A w and slopes 8 are 
obtained from the fit to (Af 2 (8)), and Eq|27l 

and the predicted signal. These deviations already sug- 
gest that our flux-selected galaxy samples are biase d w ith 
respect to the matter distribution (see Eq [T2l and lBip . 

Galaxy auto-correlation functi ons have been found t o 
closely follow a power- law (e.g. iMcCracken et al.ll2007[) . 
In order to assess the evolution of galaxy clustering with 
redshift, we fit a power- law uj{6) = A w {6/V)- & to our 
measurements, and report the best fit parameters in Ta- 
ble [H However, our estimator u>(0) is mathematically 
known only up to an additive constant known as the in- 
tegral constraint (IC), which depends on the volume and 
the m ean galaxy density of the sample (|McCracken et all 
I2007T ) . In the next section, we present how we correct for 
IC using the (Af 2 (9)} measurements. 

Galaxy-Galaxy Lensing — In Figure [6] we show the evolu- 
tion with redshift of the galaxy-galaxy lensing measure- 
ments for our 5 flux-selected galaxy samples. Overall, we 
observe our measurements to decrease with scale. To cor- 
rect for systematic effects at large scales, we subtracted 
the signal measured aroun d random foreground galaxies 
(jMandelbaum et al.l l2005h . To guid e the eye, we over - 
plot the predicted signal derived from lSmith et al.l (|2003l ) 
power spectrum with unbiased galaxies. The match is 
good at large scales, but we note a discrepancy at small 
scales. To assess the significance of this discrepancy, we 
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Fig. 6. — The mean tangential shear as a function of scale and 
redshift. The signal measured around random foreground galaxies 
has been subtracted to correct for systematic errors. To gu i de the 
eye, we overplot predictions derived from the[2mith ct al] 120031) 
power-spectrum with unbiased galaxies (solid line). We fitted 
power law functions to the data at small scales (red lines) to assess 
the significance of the discrepancy with the predictions. Error bars 
include shape noise and shot-noise (thick line) and cosmic variance 
(thin line). 
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Fig. 7. — The shear-shear correlation function for the three back- 
ground so urces catalog s B1-B 3. Solid lines are the predicted sig- 
nal with a Smith ct al. ( 2003|) power-spectrum. Error bars include 
shape noise and shot noise (thick line), and cosmic variance (thin 
line). 

compared the slopes of the data and the predictions at 
scales 9 < 3'. Taking into account the correlation be- 
tween the data points, we fitted five power-law functions 
to the data at these scales, and found their slopes to be 
always negative at 2a confidence level, hence at odd with 
the positive slopes of the p redictions. This h ighlights 
the limited accuracy of the ISmith et al.l §003) power- 
spectrum at small scales. 

Shear autocorrelation — Finally, we compute the shear 
autocorrelation functions £ + and £_ for our 3 background 
catalogs Bl, B2 and B3. They are presented in Figure [7J 
The amplitude increases between Bl and B3 mainly be- 
cause the lensing efficiency in Eq [T3] increases with red- 
shift. The contribution of the structures along the line 
of sight is likely small, because they are mostly uncorre- 
cted. 

5.2. Aperture statistics 

In this section we derive aperture statistics from the 
correlation functions. The calculations are detailed in 
Appendix [X] 
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Fig. 8. — The galaxy aperture variance. T o guide the eye, we 
overplot the predicted signal derived from a Smith ct al, II2003T ) 
power-spectrum with unbiased galaxies (solid line). As expected 
the prediction underestimates the signal at small scales. At large 
scales, the deviation in bins F4 and F5 already suggests that galax- 
ies are biased (see also Figure [TTJ. Error bars include shot- noise 
(thick line), and cosmic variance (thin line). 



Galaxy aperture variance — The measurements of 
(JV 2 (9)) for the 5 foreground bins are shown in Figure [5J 
For each redshift bin, we fit a power-law. In contrast to 
oj(9), (J\f 2 (9)) is not affected by IC, because the c ompen- 
sated filter T + which multiplies ui(9) in E q IA10I c ancel s 
out any con s tant, including IC (jSchneider et al.ll"l998[) . 
I Simon et al.l (|2007f ) derived an analytic relation between 
the parameters of a power-law fit to (AT 2 (9)) and the pa- 
rameters of a PL fit to an IC corrected correlation func- 
tion uj(6). The slope of co(9) corrected for IC is given 

by 



f(S) = o.oosia 1 



0.2769(5 



3.95 



• 0.2838a 1 



(27) 



where S is the slope of the PL fit to (M 2 {6)). The am- 
plitude of the PL fit to oj(9) is the same as the ampli- 
tude of the PL fit to (Af 2 (9)). In Tabled we report the 
IC-corrected amplitudes A w and slopes S of the galaxy 
autocorrelation function oj(9) obtained with Eq |2"T1 

We find that the amplit ude A,„ decreas e s with red- 
shift, in agreement with McCracke n et al.l (|2007D who 
also found smaller amplitudes for fainter galaxies in COS- 
MOS i.e. more likely located at higher redshift. Besides, 
we note a jump in amplitude between bins F3 and F2, i.e. 
between redshift z ~ 0.52 and z ~ 0.35. With respect 
to the slope 5, we observe a slight but not significant 
(less than la) steepening with redshift. Such a pecu- 
liar behavior could be explained by the large scale struc- 
ture identified at reds hift z ~ 0.7 already detected in 
several other studies dMassev et al.l l2007at iGuzzo et al.l 
[20071 iMeneiix et all [20091: Ide la Torre et al.l I2010D . We 
guess that more massive galaxies in this structure would 
increase the slope because they are typically more clus- 
tered than average. 

Finally, we observe an increasing deviation (more than 
5cr in bin F5) b etween the measure ments and the signal 
predicted with a lSmith et al.l (|2003| ) power spectrum and 
unbiased galaxies. This suggests that galaxies are more 
biased at high redshift, in agreement with expectations. 
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Fig. 9. — The galaxy-mass aperture covariance as a function of 
scale and redshift. To guide t he eye at larg e scale , we also show the 
predictions derived from the Smith ct al. ( 2003) power-spectrum 
with unbiased galaxies (solid line). For bin F2-B2, we found the 
change of slope of be not significant compared to a power-law fit 
(blue line). In red, (NM X ) quantifies the amount of "B-modes" 
(see text). Error bars include shape noise (thick line) and cosmic 
variance (thin line). 



For bin F4 in the last 3 bins, we find on average a bias 
of about 1.4 between the data and the predictions, and 
for bin F5 in the last 5 bins, we find a bias of about 2.6. 

Galaxy-mass aperture covariance — The measurements 
of {N{9)M &V (9)) for our 5 foreground redshift bins are 
shown in Figure [9] In agreement with the < 7t_> 
measurements, we again note that a lSmith et all ([2003D 
power-spectrum under-predicts the data points at small 
scales. We performed a x 2 -test, and found this disagree- 
ment to be significant in all bins at more than 3er (we 
included the covariances matrices in the x 2 ~test). This 
disagreement is also present for the stellar-mass selected 
galaxy samples in Figures [16] and |T7l At small scales in 
bin F2-B2, we also note a slight change of slope at 1' 
scale. Again, we investigate the significance of this fea- 
ture with a x 2 -test, but found a reduced x 2 = 1-60 which 
is not enough to reject a simple power-law model. 

An E/B mode decomposition also exits for the 
(Af(9)M ap (9)) statistics. (NM X ) quantifies the amount 
of B-modes in the measu reme nt, and is computed by re- 
placing 7 t by j x in Eq IA11I j x is the galaxy-galaxy 
clustering signal obtained with the coordinate system ro- 
tated by 45 degrees. This operation is commonly used 
to reveal B-modes contamination if a (NM X ) signal is 
detected. Using a x 2 -test and taking into account the 
correlation between the data points, we find the B-modes 
for (NM X ) to be consistent with zero in all bins Fi-Bj at 
95% confidence level, hence confirming a very low level 
of contamination in our analysis. 

Mass aperture variance — The measurements of (M 2 p (9)) 
for bins Bl, B2 and B3 are shown in Figure [TO] The solid 
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Fig. 10. — The mass-aperture variance as a function of scale and 
red shift. To guide the eye, we overplot predictions derived from 
the ISmith etaLl (2003) power spectrum. In bin F2-B2, we fitted 
a power-law (blue line) to the data and found the change of slope 
of be not significant, (in red) The "B-modes" are consistent with 
zero at all scales in all bins. Errors include shape noise and shot 
noise (thick line), and cosmic variance (thin line). 
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Fig. 11. — Evolution of bias in comoving scale and redshift for 
the flux-limited galaxy sample. Bias increases with redshift both 
because flux-selected galaxies reside in more massive halos, and be- 
cause halos of any given mass are more biased at high redshift. At 
small scales, the change of slope can be interpreted as the transi- 
tion between the two-halo and the one-halo terms in a halo model 
framework. Errors bars include shape and shot noise (thick line), 
and cosmic variance (thin line). The horizontal dashed line marks 
the val ue b = 1 for each bin. The data points are correlated (see 
Section 15.611 . 



line re presents the predicted signal with a ISmith et all 
(2003) power spectrum. We performed a % 2 -test between 
the data points and the predicted signal and found that 
the data points are consistent with the predicted signal 
at 68% confidence level for all bins Bl, B2 and B3. At 
all scales, shape noise dominates over cosmic variance. 
Errors are smaller for bin B3, because the signal is larger 
(about 10 times larger than the signal in bin Bl). 

5.3. Bias of the flux- selected sample 

5.3.1. Constant bias model 

In Figure [TT] we show the galaxy bias (6) (R) for our 
sample of flux-limited galaxies. We observe that bias 
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Fig. 12. — Same as Figure im but for stellar-mass selected galax- 
ies. Bias increases with redshift and stellar-mass. The horizontal 
dashed line marks the value 6=1 for each bin. 

varies with scale and redshift, but since the data points 
are correlated we use a x 2_ test to quantify the signifi- 
cance of these variations. First, we assume the null hy- 
pothesis that bias is scale and redshift-independent. To 
test this hypothesis, we fit a constant 6o to the 45 data 
points in the 5 redshift bins, and compute the \ 2 statis- 
tics 

X 2 = (B-6o) T C- 1 (B-6 ), (28) 

where B is a vector containing the 45 data points, and C 
is their covariance matrix. To perform the fit and find the 
best fit parameters, we use the IDL AMOEBA technique 
(jNelder fc M ead 19651), an( t repeat the process 100 times 
with different starting values. We obtain a best fit with 
X 2 = 229. Our data points are correlated, so the effective 
number of degrees of freedom (dof) to perform the % 2 test 
is less than the number of dat a points N = 45 minus 
the number of free parameters. iBretherton et all (1999) 
proposed several estimators for the number of degrees 
of freedom for correlated data. We use the following 
estimatoiP] 



dof: 



N 



E 



r 2. 



(29) 



where 



is the correlation coefficient between data 



points hi and bj, and N is the number of data points. 
Using this estimator, we find dof = 32 — 1 = 31, where 
we subtracted 1 for the parameter we fit. The reduced % 2 
is therefore % 2 /dof = 7.4, which allows us to confidently 
reject a constant bias model given our data. 

5.3.2. Redshift- dependent model 

Here we discuss a test of the redshift dependence of 
bias. To our model, we add the following redshift depen- 
dence 

9 Brcthcrton ct al. (1999) mainly studied another estimator 
dof = (y"]j On) 2 / T^.jj C?ji which is equivalent if the data points 
are centered and normally distributed with unit variance. 
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b{z) = b + b lZ . (30) 

We obtain a reduced % 2 /dof = 2.5. Although still not 
a good fit, this model is nonetheless significantly better 
than the previous redshift-independent model. 

So far, we have ruled out a constant bias with redshift 
and scale, but we did not isolate the redshift dependence 
or scale dependence yet. This is the purpose of the rest 
of this section. 

5.3.3. Scale-independent model 

Next, we test the scale-dependence only. Our null hy- 
pothesis is now that bias is scale-independent. For each 
redshift bin, we fit a constant, and sum the x 2 for each 
individual redshift bin. We obtain a total x 2 = 43 for an 
effective number of degrees of freedom dof = 32 — 5 = 27. 
This means that there are 97.4% chance that the scale- 
independent model is wrong. Here we subtracted 5 for 
the 5 constant parameters we fit. A fit of this model to 
the stellar-mass selected galaxy samples brings less strin- 
gent constraints, with x /dof = 1.49 and x /dof =1.0 
for the low and high stellar-mass selected samples respec- 
tively. 

5.4. Bias of the stellar-mass selected sample 

In Figure [121 we show the galaxy bias (b) for two 
stellar mass-selected galaxy samples as a function of 
scale and redshift. Previous studies have shown that 
stellar-mass is a good tracer of halo mass (see e.g. 
iLeauthaud et al.ll2011bl and references therein), which 
in turn parametrizes most of the bias models. In con- 
trast, flux-selected samples are more affected by selection 
effects. For instance, low surface brightness extended 
galaxies are systematically under repres ented because 
they tend to evade the magnitude cuts (jMeneux et al.l 
2009). Our stellar -mass selected samples should there- 
fore tell us more about halo properties. 

5.4.1. Halo mass 

First, we derive the average halo mass for our two sam- 
ples of galax ies. For this purpos e, we fit the bias model 
proposed bv lTinker etaLl (|2010T > 

b(u) =1 — A ^ - + Bv h + Cv c , (31) 

defined in terms of the peak height v = 5 c (z)/ D(z)a(M), 
where 6 c (z) is the critical dens ity for halo collapse 
(Weinberg & Kamionkowski 2003), D(z) is the growth 
factor, and a 2 (M) is the variance of the density fluc- 
tuations smoothed with a top-hat filter of size R = 
(3M/47rp m ) 1 / 3 . The parameters of this fitting function 
derive from fits to N-body simulations. They are A = 
1.0 + 0.24j/exp(-(4/j/) 4 ), a = O.Uy - 0.88, B = 0.183, 
b = 1.5, C = 0.019 + 0.107y + 0.19 exp(-(4/?/) 4 ) and 
c = 2.4. We define the parameter y = logioA for over- 
density A = 200 times the cosmological mean density. 

Although this bias model is not explicitly a redshift- 
dependent model, we make use of the fact that v scales 
with redshift to test the redshift dependence of our mea- 
surements. 



The iTinker et all (|2010f) model is a halo bias model, 
i.e. it was calibrated on halo clustering measured in N- 
body simulations. In contrast, we measure galaxy clus- 
tering for stellar-mass selected galaxies. Using it to infer 
halo mass implies the two following assumptions : (i) 
the dark-matter mass function is constant in the range 
of halo mass we consider (ii) there is only one galaxy 
per halo. The first assumption is valid because our 
galaxies are embedded in halos with peak height v <~ 1 
and mass log 10 (M/M Q ) < 13, a nd the halo mass func- 
tion in this regime is almost flat (|Sheth fc Tormenl 1999; 
ITinker et alJl20 "08). The second assumption holds in the 
linear regime where most of the signal comes from the 
clustering of central galaxies. 

To overcome the problem of having data points at dif- 
ferent redshifts, we fit a redshift-normalized peak height 
parameter vq, defined such that b{v{) — b(yo/ D(zi)), 
where Zj is the redshift of a data point fo;. Since this 
model is only valid in the linear regime, we only consider 
the 16 data points at scales R > 2 hr 1 Mpc. We find 
X 2 = 3.13 and % 2 = 4.93 for dof = 10.8 and dof = 11.2 
respectively for the low and high stellar-mas s galaxy sam- 
ples. Therefore, the model proposed by ITinker et al.l 
(201(J fits well our data. 

This implies that our measurements agree with an 
increase of bias with redshift. In addition, we can 
derive the halo mass of our two stellar-mass selected 
galaxy samples. Indeed, the best fit values are vq — 
0.77tg; 2 ; and u = 1.01±g;?| for the low and the 
hig h stellar-mass samples r espectively. To compare 
to ILeauthaud et al.l (|2011bf) . we compute the peak 
height estimator at redshift z = 0.37, ^ 2 =o.37 = 
^o/Z3(0.37) = 0.93 and ^ z =o.37 = 1-21 for the low and 
high stellar-mass samples, respectively. This translates 
into 9.64 < log w {M 2 oo/h- 1 M ) < 12.29, and 11.51 < 
log 10 (M 2 oo/h~ 1 M (I) ) < 12.80, respect ively. These re- 
sults ar e in very good agreement with ILeauthaud et al.l 
(|2011bf ). when considering that our two stellar-mass se- 
lected samples range between 10 9 h~ 2 M Q < M* < 
10 10 h~ 2 M Q and 10 10 fr- x M < M* < 10" fo- 2 M M . 
As a last check, we also fit the iSheth et alJ (|2001[ ) bias 
model, and found similar results. The data are not strin- 
gent enough to distinguish the two models. 

In Figure 1 1 31 we show the evolution of bias with 
redshift. The data points are averages of points at 
scales R > 2 /i _1 Mpc, and errors come from the re- 
spective total covariances matrices. The total covari- 
ance matrix is the sum of the data covariance matrix 
and the simulation covariance matrix. The shape noise, 
shot noise and cosmic variance are thus taken into ac- 
count. For c o mpar ison, we overplot the bias measured by 
IKovac et all (12 0091 (K0 9) wit h the zCOSMOS dataset, 
and bv iMarinoni et al.l (j2005f ) (M05) with the VLT VI- 
MOS Deep Survey (VVDS). The zCOSMOS dataset is 
the same as ours but with galaxies brighter than Ig,u < 
22.5. The bias measured by K09 follows our high stellar- 
mass galaxy samples, whereas the bias measured by M05 
follows our low mass samples. We attribute the differ- 
ence in bias between the two measurements to differences 
in galaxy selection functions. 

We found bias in bin F3-B3 for the high stellar-mass 
galaxy sample to be larger than expected. Although this 
data point still agrees with the bias model predictions at 
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Fig. 13. — Redshift evolution of bias averaged on scales R > 
2 h~ Mpc for the two stellar mass selected galaxy samples. Er- 
ror bars include shape noise, shot noise (thick line) and cosmic 
variance (thi n line) . Solid lines are the best fit predicted by the 
ITinker et all H20101) halo bias fitting function. At z = 0.37, the 
halo masses predicted by the model are log(Af20o/' 1 _1 Mq) = 11.58 
and log(Af20o/' l— 'Mg) = 12.36 for the low and high stellar- 
mass selected galaxy samples. Green boxes are measurements in 
VVDS for a volume-li mited (Mb < — 20 + 51og/\) galaxy sample 
IIMarin oni et al. 2005). Pink bo xes are measureme nts in zCOS- 
MOS for M B < -20 - z galaxies IIKovac et al.ir2009t) . 

95% confidence level, it shall be discussed further. In- 
deed at a similar redshift, K09 also identified a signifi- 
cant overdensit y of galaxies which c ould e xplain the large 
bias value, but iFinoguenov et al.l (|2007t) identified very- 
few massive groups. These two inconsistent observations 
are puzzling. A halo m odel alike the one proposed by 
iLeauthaud et al.l (|2011aD could help us understand bet- 
ter the peculiar properties of the galaxies at this redshift, 
in particular the fraction of satellite galaxies. 

5.4.2. Scale- dependent model 

In section I5.3[ we showed that a scale-dependent bias 
model was preferred but still not significantly better than 
a scale-independent model. Nonetheless, we observe by 
eye that bias systematically decreases at small scales, and 
the turn-down seems to start at scale R ~ 2 /i _1 Mpc, 
which could correspond to the transition between the 
one-halo and the two-halo term, alr eady identified insey- 
eral previous papers. According to Zeh avi et al.l (|2005f ). 
the behavior of bias at this transition scale is due to 
pairs of satellite g alaxies. With a halo model applied 
to simulated data, iKravtsov et al.l (|2004l ) showed that 
a significant amount of pairs wit h satellite galaxies wa s 
crucial for a smooth transition. ICoupon et al.l (120 111) , 
using the CFHT Lensing Survey, and lPeng et al.l {2011), 
with SDSS data, also found that the fraction of satel- 
lite galaxies depend on the overdensity and not on the 
halo mass or redshift. Since overdensities increase at late 
time, the number of satellite galaxies increases in low red- 
shift samples of galaxies. A change of slope in our data 



TABLE 3 

Best fit parameters for the scale and redshift dependent 
model (eq. [32l applied to the 
10 9 h- 2 M Q < M* < 10 10 /i" 2 M (low) and 
10 1() h~ 2 M & < M* < 10 11 h~ 2 M@ (high) stellar-mass 
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would therefore indicate a change of the satellite fraction 
in the galaxy samples. 

In order to measure the significance of an evolution of 
the slope with redshift, we parametrize it as a{z) = a% + 
a%z. We also introduce a turn-down scale i?TD beyond 
which bias evolves linearly. The final scale-dependent 
bias model we fit to the data is 



b(R,z) 




a 1 +a 2 z 



R < i?TD 
R > i?TD 



(32) 



w here we use t h e red shift-dependent bias model &Un( l/ ) 
of ITinker et al.l ((20100 to describe the bias in the linear 
regime. We find x 2 /dof = 0.7 and x 2 /dof = 1.0 for the 
low and high stellar-mass selected samples respectively. 
The best fit parameters are reported in Table |3] 

The difference in x 2 between the scale-independent 
(SI) and scale-dependent (SD) model Ax 2 = Xsi/dofsi — 
XsD/dofsD — 0-79 i s n °t large enough to completely rule 
out the SI bias model. 

The best fit parameters of the SD model merit some 
particular attention. First, the scale at which bias 
turns down is detected to be between 0.8 h Mpc and 
3.8 ft. Mpc at 68% CL. In contrast to what was ob- 
served in SDSS ([Johnston et al.l[2007l ). the data suggest 
i?TD to be marginally larger for less massive galaxies, but 
we would need more bins in stellar-mass to confirm this 
tendency. 

Second, the slope of bias at scales below Rtd is de- 
tected to be larger than zero, but at less than 2cr CL. 
We measured a x = 0.421°;^ and a Y = 0.63 ± 0.18 
for the low and high stellar-mass samples, respectively. 
Regarding its evolution with redshift, we found ai = 
— 017+0. 44_o. 41 for the low stellar mass galaxy sample, 
hence no significant evolution, but «2 = —0.44 ± 0.24 
for the high stellar-mass selected sample, hence a slight 
flattening of the slope at small scales at higher redshift 
bins. 

It is expected that the steepness of the slope is related 
to the occurrence of satellite galaxies in our sample, with 
more satellite galaxies producing shallower slopes. The 
coincident fact that the turn-down scale is larger and the 
slope is shallower for the low stellar-mass sample suggests 
a larger amount of satellite galaxies in our low stellar- 
mass selected sample. 

5.5. Correlation coefficients 

In Figure [T4] we show the correlation coefficient (r) 
for our flux-limited galaxy samples. Our measurements 



13 




2 

1 


2 

a i 

v 

2 



o 

+-> 

cS 

|0 

O 

o 



F5-B3 




<r> = 1 . 


3810.47 






F4-B3 


<r>= 1.3310.39 










F3-B3 














<r> = 0.58±0.28 



1 — ■-- 



.S 10 



2 

1 


2 

1 


0.1 1.0 10.0 

eff comoving scale [h _1 Mpc] 



Fig. 14. — The correlation coeffcient (r) measured in rcdshift 
bins for our flux-limited galaxy samples. A deviation from r = 1 
suggests the bias relation is stochastic, and not perfectly linear. 
Our method produces artificial deviations from r = 1, which we 
quantify with pure dark matter simulations (pink area) for which 
we know a priori that r = 1. The size of the pink area is due to 
cosmic variance. The hatched area individuates the range of scales 
where the simulations are consistent with r = 1. The averaged 
values of r are computed in these areas. 



agree with r — 1 at all scales and all redshifts. We 
nonetheless observe some trend in the data, that might 
arise from possible artifacts in the method (see Sec- 
tion IB]) . To highlight them, we applied the method to 
the simulated data, and indeed found a trend with a sig- 
nal greater than 1 at small scale and lower than 1 at 
large scale. We show the correlation coefficient and er- 
rors obtained from the simulations as pink regions. Con- 
sequently, the trends we observe in our data points can be 
explained by artifacts in the method. The pink hatched 
regions mark the scales where r is consistent with 1 in the 
simulations. To compute the correlation factor for a par- 
ticular galaxy sample, we only average the data points 
in these regions. 

We found no significant deviation from r = 1, although 
in bin Fl- Bl, we obtain r = .58 ± 0.28. IHoekstra et al.l 
(|200lD and lSimon et al.l ()2007fl found the correlation co- 
efficient (r) < 1 at 8cr, highlighting some stochasticity in 
the bias relations of their galaxy samples. In our case, 
we can only say that our flux-selected and stellar-mass 
selected samples are good tracers of the underlying mass 
distribution given the error bars (see also Figures [TBI and 
[T7]for the results with the stellar-mass samples). 

We have found no standard halo model to predict the 
scale- dependence of the cor relation coefficient r. In con- 
trast, iNevrinck et "a l. (2005) propose a galaxy-halo model 
which does. In their model, they assume halos attach to 
galaxies, in contrast to the standard halo model where 
galaxies attach to halos. By construction, assuming 
all halos have the same density profile, the correlation- 
coefficient is always r = 1. In the case halos have differ- 
ent density profiles, the correlation factor becomes scale 
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Fig. 15. — Left : The correlation matrix of estimator (b)(8) for 
data points in Figure [Til in bin Fl -Bl . Right : The correlation 
matrix of estimator (r)(8) in Figure ITU in bin Fl-Bl. Both correla- 
tion matrices include shape noise, shot noise and cosmic variance. 

dependent. At large scales, many halos are averaged 
over, resulting in a mean density profile. The matter dis- 
tribution therefore matches the galaxy distribution, and 
r = 1. At scales smaller than the minimum intergalactic 
distance, r drops below 1 because of the many different 
density profiles. Finally, assuming that the inner part of 
the density profiles are similar, r raises back to 1. 

In light of this model, we attempt to identify a dip 
in bins Fl-Bl and F 2-B2 at scale ~ 2 h ' 1 Mpc. 
IHoekstra etaH (|2002D and iSimon et all (|2007D also ob- 
tained such a dip at scale R ~ 1 h' 1 Mpc with much 
higher significance. We fit a constant r = 1 to our data 
points in the hatched areas, and measured x 2 = 3.45 
in bin Fl-Bl (dof = 2.75, > 68% CL at X 2 > 3.20), 
and x 2 = 2.20 in bin F2-B2 (dof = 3.06, > 68% CL 
at x 2 > 3.57). A constant model r = 1 therefore still 
provides a good fit to the data. 

5.6. Bias correlation matrix 

Figure [15] shows the correlation matrix of the bias 
(including cosmic variance) for the flux-limited galaxy 
sample Fl-Bl. The data points are significantly corre- 
lated between 1' and 10' for the bias parameter (b), and 
less correlated for the correlation coefficient. The large 
amount of correlation at the smallest scales 6 < 0.4' is 
due to numerical artifacts (see section [6]). 

Note as well that the large correlation between the an- 
gular bins for parameter (b) only shows up when the co- 
variance matrix derived from simulations is added to the 
covariance matrix derived from the data. Otherwise, the 
amount of correlation is weak because of the low signal- 
to-noise in the data. The signal to noise in the simulated 
data is much higher, hence the stronger signal in the cor- 
relation matrix. 

6. DISCUSSION 

6.1. Artifacts in the measurements 

It is well-known that aperture statistics, obtained 
through integration of oth er estimators, a r e very sensi- 
tive to integration limits. iKilbinger et al.l (|2006l ) found 
that the signal is underestimated by more than 10% at 
12 times the lower integration limit, i.e. 0.6' in our case. 
On the other hand, at large scales, the signal is only valid 
up to half the upper scale limit (see Eq. (9HE|) . We used 
simulations to assess the systematic effects produced by 
aperture statistics. In Figure [TBI and [T71 our simulations 
show that the correlation coefficient (r) is overestimated 
by 10% to 50% at scales 9 < 1', depending on the redshift 
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bin. For the bias parameter (b), deviations from b = 1 
occur at scales 9 < 0.7'. 

In the simulations, we also observed that the cross- 
correlation signals (M(9)M ap (9)) and (-ft) are biased low 
at large scales. Subtracting the signal measured around 
random foreground did not help in recovering r = 1 at 
large scales. We found this effect to be less important 
as the number of foreground galaxies decreases with re- 
spect to the number of background galaxies. This can be 
seen by comparing the red dashed lines in Figures IT6l and 
IT7l Although this effect is of minor consequences on our 
measurements given the size of the error bars, we tried to 
understand it. It is possible to demonstrate analytically 
that a constant additive systematic will cancel between 
sources 90 degrees separated. When foreground galaxies 
do not have sources all around them or the additive sys- 
tematic is not constant, the a verage (7^) is biased low. In 
the Sloan Digital Sky Survey, iMandelbaum et all (|2005l ) 
found that subtracting the signal measured around ran- 
dom foreground galaxies was effectively removing this 
noise. However, our simulations do not include shape 
noise, and therefore the signal we measure at large scales 
cannot be due to shape noise. 

Another source of systematic error could be related to 
the tree-code we use to compute the correlation func- 
tions. The documentation for the tree-code Athena 
mentions that choosing a too large opening angle can 
smear out the ellipticities of galaxies at large scales. In- 
deed, we find that if we increase the opening angle then 
the effect worsens, but it does not improve with values 
smaller than the one we use. Because this systematic 
error is subdominant to our measurements, we decided 
not to further attempt its correction. 

6.2. Bin mismatch in redshift 

In this section, we show that a mismatch in redshift be- 
tween peaks of the Pf(z) and the lensing efficiency curves 
g(z) can alter the bias measurements. We looked into 
this issue because this is the case in bins F4-B3 and F5- 
B3. To understand the origin of the problem, we should 
recall that in our method, the bias is the product of a 
function f% and a measurement (N 2 (9)) / (M 2 p (9)) . 

On the left panel of Figure [T8l we compute the ra- 
tios {N 2 (9)) J (M 2 p (9)} an d the inverse of the functions 
/1 for different combinations of foreground and back- 
ground sample bins. We find that between bins Fl-Bl 
and F1-B3 (Af 2 (9)) / (M 2 p (9)) decreases by 37% and 33% 
at large and small scales respectively, whereas function /1 
decreases by 78% and 64%. As a result, our estimations 
of bias change. 

We do expect the signal to decrease because (M 2 p {9)) 
increases with redshift, but we do not expect the mea- 
sured and the predicted signal to change by different 
amounts. According to Eq 1191 the bias measurements 
should only depend on the foreground redshift distribu- 
tion, and not on the choice of the background galaxy 
sample and associated lensing efficiency. This is of course 
assuming that (M 2 p (9)) and the matter power spectrum 
grows linearly, and the lensing-effective matter redshift 
distribution is well approximated by the lensing effi- 
ciency. For very wide surveys where many line of sights 
are averaged over, the lensing efficiency might be a good 
approximation of the effective matter distribution, but 



in COSMOS especially at low redshift, this might not 
be the case. On the other hand at small scales in the 
non-linear regime, the power spectrum might not grow 
linearly. In Figure HH smaller scales are more affected 
by bin mismatch. 

Given the limited size of the COSMOS held, we 
therefore posit that cosmic variance or shape measure- 
ments could yield these discrepancies because the vol- 
umes probed by (M 2 p (9)) are too small. This type of 
discrepancies warn us that for future surveys, we should 
try to match as closely as possible the mean redshift of 
the foreground galaxies and the peak of the lensing effi- 
ciency, in order to limit the effect of cosmic variance. 

7. CONCLUSION 

The strength of the COSMOS survey is the exceptional 
quality of the ACS imaging, and the 30-band photome- 
try. Thanks to the latter, a photometric redshift can be 
derived for more t han 600,000 galaxies at Isuw < 26.0 
(|Ilbert et al.l[2009h . 

In this paper we made use of these two properties to 
investigate the evolution of stochastic bias with scale 
and redshift in COSMOS. We partitioned our foreground 
galaxy catalogs in 3 categories : one flux-selected cat- 
alog (Isuw < 26.5 and K s < 24), and two stellar- 
mass selected catalogs. To estimate the bias parameter 
b and the correlation coefficient r, we appli ed a tech- 
nique based on aperture sta tistics described in (Schneider 
Il998t iHoekstra et al.l 120021 ) . We used simulated lensing 
catalogs derived from N-body simulations to quantify 
the amount of error due to cosmic variance, and to test 
the method against numerical artifacts. We found that 
weak lensing shape noise dominates the error budget, and 
that we are affected by some numerical artifacts at small 
scales (9 < 0.6'), and at large scales for the correlation 
coefficient r. 

Then, we used a x 2_ test to assess the significance of 
redshift evolution and scale-dependence of bias. We 
found a significant evolution of bias with redshift, 
but could not confidently rule out a scale-independent 
model. We f ound that both bias models p roposed by 
iTinker et all (pOlOt ) and iSheth et al.1 (pOOlf l provide a 
good fit to the redshift evolution of our measurements 
at scales R > 2 /i _1 Mpc. We used the result of the fit to 
derive the halo mass of our stellar-mass selected galaxy 
samples. Finally, we observed that bias starts to decline 
below a scale of about 2 ft, _1 Mpc. We proposed a bias 
model to describe this scale-dependence and obtained a 
very good fit with x 2 /dof < 1 for our two stellar-mass se- 
lected galaxy samples. We measured a 3a significance of 
the turn-down scale Rtd for the high stellar mass sam- 
ple, but could not draw conclusions about the evolution 
of the slope with redshift because of low signal to noise. 

Finally, we observ ed that bias in bin F 3-B3 is 2a off the 
expectations with a lTinker et al.l (|2010D bias model. We 
attribute this to the presence of large over-densitie s in the 
CO SMOS field al r eady id entified in IKovac et al.l ff2009h 
and lMassey et alJ (|2007al ). With respect to bias stochas- 
ticity parameter (r) , our measurements are all consistent 
with (r) =1. The linear bias model is therefore a good 
fit to our data, albeit with somewhat large errors. 

The COSMOS field is still too small to produce deci- 
sive conclusions about the evolution of bias, especially 
to characterize its scale-dependence, and to measure its 
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stochasticity in stellar mass-selected galaxy samples. On 
average, we measured signal to noise ratios bja^ ~ 3 and 
r/a r ^0.5. The most important source of error is shape 
noise in the background galaxies. iSchneider et al.1 ()1998[ ) 
have shown that the error on (Mf p (0)) decreases as 1/Nf,. 
In order to increase the signal to noise by a factor of 10, 
we would need a field 10 times larger. 

With t heir RCS and VIRMO S-DESCART survey of 
50.5 dee 2 . iHockstr a et al.l ()2002f ) obtained signal to noise 
ratios about 8 on both bias parameters b and r. However, 
they did not perform simulations to disentangle system- 
atic and statistical errors as we did in this paper. They 
also did not have redshift or stellar-mass estimates for 
their galaxies. Future ground-based wide field surveys 
such as KIDS and VIKING, DES, Pan-STARSS, LSST 
and HyperSuprime-Cam on Subaru should provide good 
photometric rcdshifts and stellar-mass to perform a more 
conclusive analysis of the scale- and redshift-dependence 
of bias. 

The main challenge for weak-lensing tomography with 
the future surveys (this bias analysis being one applica- 
tion) is to measure the shape of background galaxies up 
to redshift z = 4. This is essential to measure the red- 
shift evolution of galaxy properties and bias up to red- 
shift z ~ 1. In their current implementation, the space- 
based missions Euclid and WFIRST propose a near- 
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APPENDIX 
PRACTICAL ESTIMATORS 

In practice, aperture statistics estimators are obtained from two point correlation functions (Schneid er et al] |2002) . 
The correlators used in this work are : (i) the angular correlation of the foreground galaxy positions, w(0), (ii) the mean 
tangential shear about foreground galaxies, (7t)(0), (hi) the sh ear-shear correlation functions £,±(0) as determined from 
the ellipticities of the background galaxies (jSimon et~aT1[2007l) : 

The angular correlation of the foreground galaxy positions is computed with the lLandy fc Szalayl (|1993[ ) estimator 
defined as 

The mean tangential shear around foreground galaxies is computed in angular scales, and is function of the shear 
components 71 and 72 

(7t)(0) = {w^ i ' j) ^ i ' j) (O)) (A2) 

( 7 W cos2<^ + 7 W mu2o" ) A^) 

where <f>( % 'i' — 0W — #0) is the polar angle between foreground galaxy i and background galaxy j. The angle brackets 
stand for weighted mean over indices i and j and weig hts W W. A( i ^')(6') is a binning function 

\ otherwise 

The shear-shear correlation functions £±(0) are obtained from the correlation of the shear components 71 and 72 : 

e + (0) = / w « U7 «) 7 W 7 tj)^ + (w^w^^^) A^(0) (A4) 

£_(0) = (Jw W w ^hi Videos 4<^ 
- ( w^w^hohi j) cos 46^) 



w 



w W w 0) 7 W 7 W sin40( 4J ^) A^(0) 



Aperture statistics estimators can be expressed in terms of these correlators (jSchneider et al.ll2002t [Hockstra et al.1 
2002). First, the E- and B-mode of the aperture mass variance (M| L(0)) are obtained from the shear-shear correlation 
functions £±(0) 

<Km-ir^(t>wJi) + t-w-(i)). (A6) 



where 



1 f 2f> M 1? 



(Mlm = - — ( e + WT+ ( 7T ) ( 77 ) ) • (A7) 



T + (X): 



6(2 - I5x 2 ) 



arcsin(a;/2) 



100tt 
+132x 6 - 9x 8 ) 



(120 + 2320x 2 - 754 4 



(A8) 
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r -«- IK 1 -£)'"■ < A9 > 

T±(x) vanish for x > 2. The B-mode aperture mass (Mj_(0)) provides a quantitative estimate of the systematics, 
since gravitational lensing only produces E-modcs. 
Second, the aperture count variance (Af 2 (9)) is related to u(8) 

<^ 2 W> = ^ (£) ■ (A10) 

Finally, the galaxy-mass cross correlation aperture (Af(0)M ap (8)) is obtained from the mean tangential shear (jt)(9) 

29 dtf d „ (-a 



(Af(9)M ap (9)) = J —{ lt )(6)F [ - ) (All) 



where the function F{x) has no analytic expression and also vanishes for x > 2 



f 2 c\t 

F(x) = 576 / -J 2 (xt) [J A {t)f . (A12) 
CORRELATORS 

We have seen in Eg 151 - ITD1 that the aperture statistics {M 2 {8)), (M 2 p (6)} and (Af(8)M ap (6)} derive from their 
respective power-spectra multiplied by the same filter function I (£9) = Ji(x)/x 2 . This filter is convenient because 
it is very narrow in Fourier-space, and almost acts as a 6d(x) function. When integrated over the power-spectra, it 
results in very limited correlation between the Fourier modes. The aperture statistics are close approximations of the 
power-spectra but in real space. There is almost a bijective relation between angular scales, and Fourier modes. As a 
result, it is possible to combine the aperture statistics to derive the bias parameters as it is done in Fourier space. 

In practice, aperture statistics derive from the correlators u>( 9), ( r )ti{9) and£+(fl) . Correlators also derive from the 
power-spectra, but filtered by different filters, as shown below (Hoekstra et al. 2002) 



u>(6) = (5n(0)6n(6)) 
1 

~ 2^ 



d£ I P n (£)J o (£0) , (Bl) 



(i 



d£ I P Kn (£)J 2 {£9) , (B2) 



e±(0) = <7t(O)7 t (0))±<7x(O)7x(0)) 

= h I" M * p *(*)JoAM) > (B3) 



Therefore, taking their ratio to derive the bias parameters is unreliable. Note as well that some filters alike the Besel 
function Jo(x) act as low-pass filters. Many Fourier modes are mixed together into one angular scale. Consequently, 
the correlation between the different angular scales is large. 
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Fig. 16. — Overview of the most important estimators involved in the measurement of bias and stochasticity for the stellar mass selected 
galaxy sam pie in the ran ge 10 9 < M* < 10 10 h~ 2 M©. The black curve in the two upper panels is the theoretical prediction derived 
from a Si:::ll: v. ;,]. (2003) power-spectrum, and it is meant to guide the eye at large scales. The discrepancies at small scales are expected 
(see text). The red dashed curves are the signals obtained in the simulations, which highlight numerical artifacts. Data points must be 
interpreted with respect to these curves rather than b = 1 and r = 1. 
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Fig. 17. — Same as Figure H"6l but for galaxies in the stellar mass range 10 10 < M, < 10 1 
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Fig. 18. — The effects of a mismatch between foreground (Fl, F2, F3) and background (Bl, B2, B3) rcdshift bins for the stellar mass 
selected galaxy sample 10 9 < M» < 10 lc> . Left panel — The bias parameter (6) not calibrated by function fl. The thick line is the 
calibration signal l//i computed assuming 6=1. Right panel — The bias parameter (6) calibrated by function fl. 



